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ABSTRACT 



We calculate the orbital evolution of binary systems where the primary star 
is an evolved red giant branch (RGB) star, while the secondary star is a low- 



Csl ■ mass main sequence (MS) star or a brown dwarf. The evolution starts with 

. a tidal interaction causes the secondary to spiral-in. Than either a common 

1-^ ! envelope (CE) is formed in a very short time, or alternatively the system reaches 

\0 '. synchronization and the spiraling-in process substantially slows down. Some of 

r— |! the latter systems later enter a CE phase. We find that for a large range of system 
parameters, binary systems reach stable synchronized orbits before the onset of 

• ! a CE phase. Such stable synchronized orbits allow the RGB star to lose mass 

p ^! prior to the onset of the CE phase. Even after the secondary enters the giant 

Q ! envelope, the rotational velocity is high enough to cause an enhanced mass-loss 

^ ! rate. Our results imply that it is crucial to include the pre-CE evolution when 

c3 . studying the outcome of the CE phase. We find that many more systems survive 

! the CE phase than would be the case if these preceding spin-up and mass-loss 

^ ! phases had not been taken into account. Although we have made the calculations 

CN ! for RGB stars, the results have implications for other evolved stars that interact 

^ ! with close companions. 

o : 
d ■ 

0\ '. 1. Introduction 

O . 

' Extreme horizontal branch (EHB) stars are post red giant branch (RGB) stars that burn 

^ ■ hehum in their core and have a very low-mass envelope. Spectroscopicahy they are classified as 

■ sdB (or sdO) stars, although some sdB and sdO stars might be post asymptotic giant branch (AGB) 

stars. In what follows we will refer by sdB stars only to EHB stars, and include under the name 
sdB stars also the sdO stars. Many of the sdB stars in the field (i.e., not in globular clusters) are 
in binary systems (Maxted et al. 2001a,b), some with a close and very faint M type main sequence 
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(MS) companion, e.g., HS 0705 + 6700 (Drechsel et al 2001), PG 1336-018 (Kilkenny et al. 1998), 
and HW Vir (Menzies & Marang 1986; Wood et al. 1993). 

The orbital and physical parameters of the three systems HS 0705 + 6700, PG 1336 — 018, 
and HW Vir, are very similar (Drechsel et al 2001). In HS 0705 + 6700, for example, the sdB 
and secondary stellar masses and radii are Mi = O.483M0, M2 = 0.134Mq, Ri = O.23Oi20, 
and R2 = O.186i?0, respectively, while the orbital separation and period are O.81i?0, and 2'^18™, 
respectively. It turns out that the secondary star almost fills its roche lobe (Drechsel et al 2001). 
Typically these systems have q = M2/M1 < 0.3. Detection is biased toward close and more massive 
companions. This suggests that many more systems exist with even lower values of q (Davis et al. 
2009). Geier et al. (2009), for example, announced recently the detection of a 8 — 23Mj, where Mj 
is Jupiter mass, substellar companion orbiting an sdB star at an orbital separation of 5 — 6.1i?0. 

These systems must have gone through a common envelope (CE) phase (e.g.,Livio & Soker 
1988; Drechsel et al. 2001; Han et al. 2007 and Han 2008 for recent papers with more references). 
Different aspects regarding the formation of sdB stars and the evolution of RGB stars through 
a CE phase have been discussed in the literature for more than three decades (e.g., Pustynski & 
Pustylnik 2006; Soker 2006; Nordhaus & Blackman 2006; Nordhaus et al 2007; Livio & Soker 2002; 
Dewi & Tauris 2000; D'Cruz et al 1996; Alexander 1976; Paczynski 1971; Han et al. 2002, 2003, 
2007; Han 2008; Zuckerman et al 2008; Melis et al 2009). However, some questions remain open. 
(1) Does the companion transfer part, or all, of the mass it accreted during the CE phase back 
to the primary at the end of the CE phase? (2) What is the role played by very low-mass MS 
companions, considering that even substellar objects (brown dwarfs and planets) can form sdB 
stars, as suggested by one of us in a series of papers (Soker 1998, Soker &: Harpaz 2000, 2007), and 
confirmed observationally by Geier et al. (2009)? 

In the present paper we examine how tidal interaction before the formation of the CE might 
increase the likelihood of a low-mass MS companion to survive the CE phase. In section [2] we 
explain the basic mechanism and present simple solutions ignoring the evolution of the primary 
star along the RGB (but mass- loss is included) . In section [3] we also consider the increase in the 
stellar core mass and stellar radius. We summarize our main results in section 21 

2. Basic characteristics of the evolution 

The onset of the CE, and the ejection of the envelope before and after the onset of the CE, 
have received a great deal of attention over the years (e.g., Soker et al 1984; Tauris & Dewi 2001; 
Han et al. 2002; De Marco et al. 2003; Soker 2004; Hu et al. 2007; Nelemans & Tout 2005; Beer 
et al. 2007; Politano & Weiler 2007; Webbink 2008; Davis et al. 2009). Our goal is to study in 
more detail the evolution of binary systems in a stage prior to the onset of the CE phase, and in 
particular systems that have reached synchronization; the synchronization is between the orbital 
period and the primary rotation period. 
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We start our calculation when tidal interaction becomes important. For the binary systems 
we study, where the primary is an RGB star and the secondary is a low-mass main sequence (MS) 
star or a brown dwarf, tidal interaction becomes important when the giant swells to a radius of 
Rg ~ 0.2a, where a is the orbital separation (Soker 1998). The primary radius increases along the 
RGB as the core mass increases. The primary and the secondary initial masses range in the binary 
systems we study are O.SMq < Mi < 2.2Mq, and O.OISMq < M2 < 0.2Mq, respectively. 

When a binary system starts its evolution it is not synchronized, and therefore tidal interaction 
will lead to a fast spiraling-in process, i.e., the orbital separation decreases. The binary system 
can then either reach a synchronization or stays asynchronous. In the asynchronous systems tidal 
interaction will lead to a fast spiraling-in process all the way to a rapid formation of a CE. In 
systems where synchronization is reached the orbital evolution substantially slows down (unless 
a Darwin instability occurs as soon as synchronization is achieved; see below), as its rate is now 
controlled by mass-loss and the evolution of the primary star; the orbital separation can decrease 
or increase. 

When the orbital separation increases faster than the RGB stellar radius, the evolution is 
terminated when all the primary envelope is lost, or a core helium flash occurs; no CE phase 
occurs. If the orbital separation decreases relative to the RGB stellar radius, the evolution ends 
when at least one of the following conditions occur: 

(1) When the mass of the primary envelope has been totally lost. No CE phase occurs. 

(2) When the core experiences a core helium flash. Here we take this to occur when the core mass 
reaches O.48M0. No CE phase occurs. 

(3) Darwin Instability (see below) is reached, the secondary falls into the primary envelope. 

(4) The RGB stelar radius increases to the point where it engulfs the secondary, namely, Rg > a. 
If conditions (3) or (4) are met, the system enters a CE phase after some fraction, even a large 
fraction, of the primary envelope has been lost. 

The different types of behavior are illustrated in Fig. [T] This figure includes the evolution of 
the primary star, and the detail of the calculations will be presented in section [3J Here we are only 
interested in the general outcome. The orbital separation (in units of Rq) versus time is presented 
for several cases, as indicated in the figure. When the secondary mass is too low (exact value 
depends on the RGB stellar properties) the secondary star does not have enough orbital angular 
momentum to bring the primary rotation in synchronism with the orbital period. Tidal interaction 
rapidly cause the secondary to spiral-in and form a CE. We do not calculate the tidal time scale, 
as it is very short. 

When the secondary mass is larger than some critical value for the same RGB star, a synchro- 
nization is achieved (in a short time) before the orbital separation decreases to the RGB radius. (As 
we do not treat the initial process of reaching synchronization, the calculation seems to start from 
a radius of smaller that the primordial orbital separation.) Now the system evolves on a mass- loss 
time scale. The primary star loses mass, an effect that tends to increase orbital separation. How- 
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Fig. 1. — Orbital separation, in solar radii for different binaries as indicated, vs. time. The 
termination of the evolution is indicated. NoSync: No stable synchronized orbit is reached, 
and a rapid CE phase occurs; Menv : the RGB has lost its envelope; Rg = a : the RGB 
stellar radius increased to the point where it swallowed the secondary; Darwin: the Darwin 
instability caused the formation of a common envelope; He flash: the core mass reached a 
mass of O.48M0 and a He core flash was assumed. The initial RGB stellar mass Mg and 
secondary mass M2 for each case are indicated in the legend. In all cases the mass-loss 
parameter in equation ([8]) is r^/j = 3 x 10~^^, the initial core mass is Mc(0) = OAMq, the 
orbital separation prior to tidal interaction is = 5Rg{0) = 4O5i?0, and the calculations 
include the evolution of the primary star as explained in section [31 
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ever, the wind carries angular momentum. To maintain synchronization the secondary transfers, 
via tidal forces, orbital angular momentum to the envelope, and the orbit shrinks. Tidal effects 
might overcomes the mass-loss effect. Eventually, a point is reached where the secondary enters 
the envelope, either due to the Darwin instability, or by the swelling RGB envelope that increases 
to Rg = a. By that time the envelope mass is lower than its initial value. Calculations that do not 
consider the possibility of pre-CE mass-loss might attribute the entire mass-loss to the CE phase. 
Such calculations overestimate the efficiency by which the secondary expels the envelope during 
the CE phase. This overestimate can lead to a search for an extra energy source that is actually 
not required. 

When the secondary mass is larger even (for the same RGB star), the secondary orbit shrinks 
by a small amount before it brings the primary to synchronization. As the orbital separation is 
large, the effect of mass-loss overcomes the effect of angular momentum transfer by tidal interaction 
and the orbital separation increases. 

To demonstrate the basic characteristics of the evolution, we start by neglecting the evolution 
of the RGB star. Namely, we assume the core mass and RGB stellar radius to be constant, as in 
Soker (2002). As a single star reaches this radius its rotation angular velocity oj will be very small; 
practically zero. When synchronization is achieved, it is done on a short time, and we neglect 
angular momentum loss in the wind, such that the source of most of the RGB angular momentum 
Jenv = lenv^, IS the change in the orbital angular momentum AJorb, where /gnv is the envelope's 
moment of inertia. Namely, 

-^env'^syn — At/grb) (1) 

where Wsyn is the angular velocity of the primary envelope when synchronization is achieved for the 
first time, and 

AJorb = KGMf' (a[]-^ - flO/J , (2) 

where oq is the orbital separation prior to synchronization (when tidal interaction becomes impor- 
tant), and is set here to qq = 5Rg{0), and Ogyn is the orbital separation when synchronization is 
achieved for the first time. As usual, M = Mg + M2 and /x = MgM2/M where Mg = Mi is the 
mass of the giant. Synchronization implies the equality 

UJ,yn = [GMgia.ynr'^f^ (S) 

for M2 < Mp. 

Due to its small size, the core's angular momentum is negligible. The angular momentum 

carried by the wind must be considered after synchronization is achieved because the evolution 
slows down, and the primary's mass-loss rate on the RGB is significant over a long time. Actually, 
under the assumption that the primarys core does not evolve, the evolution is controlled by mass- 
loss. Conservation of angular momentum determines the evolution 



■-^orb "I" "-^env "I" "Avind — 0. 



(4) 
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Using the evaluation of equation Q as given by Soker (2002), we solved for the orbital evolution 
as function of mass-loss for an RGB star with constant core mass of Mc = OAMq and radius of 
Rg{0) = 81Rq. Tidal interaction is assumed to start when the stellar radius grows to Rg = 0.2oo, 
hence here we set ao = 5Rg{0) = A05Rq, and a constant core mass of Mc = OAMq (for details see 
Soker 2002). The RGB radius is fixed at Rg{0) = 81Rq. When the evolution of the primary is 
not considered, the orbital separation is a function of the envelope mass. For that case there is an 
analytical solution to determined if a stable synchronization orbit is achieved. 

In some cases when synchronization is achieved, the system is already unstable to the Darwin 
instability. This instability sets in when a small decrease in the orbital radius does not transfer 
enough angular momentum to the spinning primary to maintain a synchronization. The condition 
for stability is lorb > 3/env) where lorb = MgM2/{Mg + M2)a^ is the orbital moment of inertia, 
and /env = V^envR'^ for the envelope's moment of inertia. For RGB stars rj ~ 0.2. Taking again 
M2 -C Mg, the instability sets in when 

Maa^ = Sr^MenviJg. (5) 

Equations (H]), ([2]), ([SD and Q can be combined to yield the condition that a stable synchro- 
nized orbit is achieved, that for M2 ^ Mg can be approximated as 

Taking oq = 5Rg as before, rj = 0.2, and Mgnv = Mg — 0.4, for a core mass of Mc = OAMq (masses 
in solar units if not indicated otherwise), in equation ([6]) gives 

M2 > 0.076Menv = 0.076(Mg - OAMq). (7) 

This represents the straight line in the Mg — M2 plane, drawn as straight blue line in Fig. [21 

In Fig. [2] we present by a wide strip the region in the Mg — M2 plane where the systems 
reach synchronization before the onset of a CE. Above this strip the companion spiral in into the 
envelope before a synchronization is achieved. Below that strip, the secondary avoids the CE. As 
found already by Soker (2002), this strip is very narrow. The line according to equation ^ is also 
drawn in the figure, just on the upper-left boarder of the strip. To better explore the possibility of 
significant mass-loss prior to the onset of a common envelope, we turn to include the evolution of 
the primary. 

3. Calculations with an evolving primary star 

3.1. Basic equations 

When the evolution of the primary is not considered, the orbital separation is a function solely 
of the envelope mass, as the core mass and radius are constant. In that case the dependence on 
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M2{Mq) 

Fig. 2. — Mg vs. M2 for binaries with a constant RGB core mass M^, and hence constant 
radius Rg (RGB evolution is not considered). The upper-left trapeze represents binaries that 
did not achieve synchronization, but rather formed an immediate CE. The lower-right trapeze 
represents binaries that achieved synchronization, but never formed a CE. In these cases the 
envelope was depleted, mainly due to the RGB wind. The thick strip in the M2/Mg plane 
is for binaries that achieved synchronization before they entered a late CE. The straight 
line represents equation ([7]). In all cases here the mass-loss parameter in equation ([8j) is 
?7h = 3 X 10~^^, the core mass is Mc = OAMq, and the orbital separation prior to onset of 
tidal interaction (primordial separation) is qq = 5Rg{0) = 405/?©. 
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time comes from the mass-loss rate. In this section the increase in core mass and its influence on the 
primary RGB star are taken into consideration. The mass-loss from the envelope is calculated by 
using the Reimers equation (Rcimcrs 1975; sec, e.g., Catclan 2009, Schoder & Cuntz 2007, Suzuki 
2007, and Meszaros et al. 2009 for recent usages and discussions of that formula) 

Mg = ijrLqRq/Mq (8) 

Although this formula is commonly used to estimate the mass- loss from the giant, most researchers 
agree that it does not cover all the physical mechanism and characteristics. For example, we know 
from the different distributions of HB stars in different globular clusters that other effects, such 
as mctallicity, influence the wind ejection process (e.g., D'Cruz et al 1996, 2000; Whitney et al. 
1998; Catelan 2000). The value of r/^ has been discussed at length in the literature. Since it was 
shown that the Reimers formula underestimates the wind ejected from some giants (Schroder & 
Kuntz 2005, Suzuki 2007), we will explore the results for two values: rjR = 3x 10~^^ (normal), and 
riR = ^x 10-'^ (high). 

In this section we consider the nuclear evolution of the star. Nuclear burning increases the 
core mass, implying increase in luminosity and RGB radius. For the growth rate of the RGB core 
we take an expression given by Iben Sz Tutukov (1984) and Nelemans & Tauris (1998), based on 
results of Mengel et al (1979) 

M, = 10-5-36 /^y '^^y^-i_ (9^ 



The RGB stellar radius is calculated according to the relation proposed by Iben &; Tutkov (1984) 

The rate of angular momentum carried by the wind (defined positively) in a synchronized 
system (primarys rotation synchronized with orbital motion) is given by (Soker 2002) 



J. 



wind 



^2 y o , 



a;(-M,). (11) 



The first term is the angular momentum due to the orbital motion of the primary around the 
system's center of mass, while the second term is due to the rotation (spin) of the primary. Here 
—Mg = Myj > is the mass-loss rate in the wind. We neglect accretion by the secondary. The 
rate of change in the orbital angular momentum Jorb = //(GMa)^/^, where n = MgM2/M is the 
reduced mass, and M = Mg + M2 is the total mass of the binary system, is given by 

^ = #-0-577^+0.51 (12) 
Jorb Mg Mg + M2 a ^ ' 
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The rate of change of the envelope angular momentum Jgnv = Icnv^ is Jcnv = Icnv^ + "^-^env 
Substituting /env = il^envR'^, the synchronization condition uj = ^ G{Mg + M2) / , and assuming 
that r] does not vary along the RGB, gives after some rearrangement 



J env ^cmi . ri^g in r 



M„ + M2 ■ a ■ 



(13) 



Substituting Eqs. ([II])- ([l3|) into Eq. (gD gives (M = + M2): 



1 - 3r/Me 



1-2/3 



M„M2 V « 



-2ry 



MenvM 



MgM2 V a 



Rg 



MenvMg 
M„M2 



fig M 

R'„M„M2 



R„ 



(14) 



3.2. Solution procedure 

The solution was obtained under the following assumptions and by performing the following 
steps. 

(1) We chose as initial conditions: Core mass Mc- This fixed the RGB radius by equation (jlOp . 
RGB stellar mass Mg] companion mass M2; orbital separation a with a circular orbit; coefficient 
ijfi of the mass-loss formula ([8|). For the RGB envelope's moment of inertia we take ry = 0.2, and 
for the angular momentum loss we take /3 = 2/3 in equation (|lip . as appropriate for a uniform 
(spherical) mass- loss geometry from the surface of the RGB star. 

(2) The tidal interaction depends very strongly on the ratio of the stellar radius to orbital separation 
(Rg/a). We therefore assume, as in section [2l that tidal interaction is negligible, until the RGB 
stellar radius is large enough, when a strong tidal interaction occurs. We take this ratio to be 
Rg = 0.2a. 

(3) As in section [21 the onset of the tidal interaction results in one of three possibilities, (i) The 
orbital angular momentum is too low, and no synchronization is achieved. The companion spirals- 
in to the RGB envelope, (m) Synchronization between the orbital motion and the RGB rotation is 
reached, but the system at that stage is unstable to the Darwin instability (see eq. [5]). Here again, 
the companion spirals-in within a very short time to form a CE. (in) Synchronization is achieved 
and the system is stable. This evolutionary rout is further evolved, as described below. 

(4) After synchronization is achieved in a stable system, we start to evolve the primary star. The 
growth of the core mass is according to the rate of hydrogen nuclear burning that forms helium, as 
given in equation ([9]). This mass is removed from the envelope. It changes the moment of inertia 
of the envelope, but not its angular momentum. 

(5) Simultaneously with the growth of the core we consider the wind, by removing mass and angular 
momentum from the RGB star, according to equations ([8]) and (jlip . respectively. 

(6) As the wind caries angular momentum, and the growth of the core changes the RGB moment 
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of inertia, in each time steps we solve equation (|14p to determine a new orbital separation a. 
(7) In each time step we check whether one of the following occurs, (i) The RGB envelope has been 
completely depleted. We register the orbital separation, (ii) The RGB core has reached a mass of 
Mc = 0.48Mq, when a core helium flash is assumed to occur. We register the orbital separation. 
(iii) The Darwin instability occurs. In this case a late CE is formed. We register the RGB envelope 
rotation velocity, (iv) The increasing radius of the RGB reached the orbital separation Rg = a. In 
this case a late CE occurs. 



3.3. Results 

3.3.1. Evolution of orbital separation and envelope mass 

We solved for the evolution according to the steps outlined in section 13.21 In Figs. [3] and H] 
we present the outcome of the evolution for our two mass-loss rates: ?/_r = 3 x 10~^^ (normal) and 
r]R = 3 X 10~^^ (high), respectively, in equation In these sets of calculations the initial core 
mass is Mc(0) = 0.4Mq, such that the RGB radius is Rgo = Sli?©. As strong tidal interaction 
starts when a = 5Rg by our standard parameterizations, so that the initial binary separation is 
a = mRQ. 

For the normal mass-loss rate (Fig. [3]) the envelope mass decreases slowly, and so does its 
moment of inertia. As the RGB radius increases, the Darwin instability sets in and the companion 
spirals-in to the RGB envelope to form a CE. The CE outcome covers a large area in the M2 — Mg 
plane. It is interesting to compare this case to the case where the RGB evolution is not considered 
(Fig. [21 Soker 2002). When the increase in core mass and radius are not considered, only a 
small parameters-area in the M2 — Mg plane leads to the formation of a CE phase after a stable 
synchronization has been achieved. When the increase in radius is taken into account, this area 
is much larger. For a small area in the M2 — Mg plane (lower-right) the core reaches a mass of 
Mc = O.48M0, for which a core helium flash is assumed to occur (turning the star to a horizontal 
branch star), and the binary system does not enter a CE phase. Only much later during the AGB 
phase the system might enter a CE phase. 

For the high mass-loss rate (Fig. the envelope mass decreases more rapidly, and so does 
its moment of inertia. This makes the Darwing instability less likely. For that, in some cases the 
RGB radius reaches the orbital separation before a Darwin instability occurs. Also, the higher 
mass-loss rate leads to the depletion of the envelope before a Darwin instability occurs in some 
cases. In addition, mass loss acts to increase the orbital separation. Therefore, for a larger range 
of parameters-area in the M2 — Mg plane (lower-right corner) the system does not enter a common 
envelope at all. Still, binaries that located in a large area in the M2 — Mg plane do enter a CE 
phase after a stable synchronization has been achieved. For these systems the RGB envelope mass 
decreases, in some cases substantially, before a CE phase sets in. 
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Normal mass loss rate 




Fig. 3. — Mg vs. M2 representing binary systems that reached or did not reach synchro- 
nization. Our normal mass-loss rate of 77/^ = 3 x 10~^^ was used. The initial core mass is 
Mc(0) = OAMq, and the initial (primordial) orbital separation prior to tidal interaction is 
Qq = 5Rg{0) = 4O5i?0. Calculation is terminated as indicated, when one of the following 
occurs. (1) The Darwin instability brings the system to a CE phase (marked 'Darwin Insta- 
bility'); (2) The RGB stellar radius exceeds the orbital separation {Rg = a; does not happen 
in this case); (3) A total depletion of the RGB envelope occurs (Menv; does not occur in this 
case due to a relatively low-mass loss rate as indicated by rj^). (4) The core mass reaches 
0.48Mq and assumed to go through a core Helium flash ('Helium flash'). 
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High mass loss rate 




M2{Mq) 



Fig. 4. — Like Fig. |3l but for a high mass-loss rate with rjn = 3 ■ 10 
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In Fig. Owe present the evolution of the scaled orbital separation, a/ao, with time, for the 
two values of the mass-loss rate, and for several RGB initial masses. In all runs = 405/^0 and 
M2 = 0.2Mq, and all cases that are shown reach a stable synchronized orbit after the initial tidal 
interaction. The following characteristic of the evolution should be noted. 

(1) The values of a/ao start with values of a/ao < because we assume that the initial strong 
tidal interaction brings the system to synchronization in a short time. 

(2) Since we start the calculation when the RGB is highly evolved (a core mass of Mc(0) = OAMq 
and a radius of Rgo = 81Rq), the evolution to the end of the RGB is relatively short, < 4 x 10^ yr. 
If the initial orbital separation is smaller, then strong tidal interaction occurs earlier, when the 
evolution is slower, and the evolution time will be longer. 

(3) Synchronization implies that the RGB star rotates at about ^0.1 times it break-up speed. 
Namely, ^0.1. With its strong envelope convection, the synchronized RGB star is likely to have 
strong magnetic activity. The rotation and magnetic activity are likely to substantially enhance 
the mass-loss rate. For that, we take the high mass-loss rate cases to represent our expectation for 
the mass-loss rate. 

(4) As expected, a higher mass-loss rate reduces the chances for the formation of a CE. Still, by 
Fig. m for a large parameters-area a CE will occur after a substantial mass of the envelope has 
been lost. The decreases of the envelope mass with time until the onset of a CE phase is presented 
in Fig. [6l The envelope mass is scaled with its value at the beginning of our calculation (when the 
system reaches synchronization for the first time). 



We saw that many systems form a CE phase only after the envelope mass has been substantially 
reduced. Another important effect in all the evolutionary routes considered here is the spinning-up 
of the RGB-star envelope before the CE phase starts even. With their strong convection, RGB 
stars that rotate even at a moderate rate can have a relatively very high mass-loss rate. 

The envelope angular velocity LVg should be compared with the break up velocity (the velocity 
of a test body performing a circular Keplerian orbit on the stellar equator) 



The important quantity here is the ratio ilg = ujg/ujB- 

The orbital velocity of the RGB envelope ujg changes during the evolution. We will consider 
four points in the possible evolutionary routes: Vtgo : for systems that did not reach a stable 
synchronized orbit, this is the angular velocity of the envelope when the secondary just enters the 
envelope; Vlgi : when stable synchronization is first achieved; : when a system with a stable 
synchronized orbit encounters the Darwin instability and loses synchronization; ; the primary 



3.3.2. Envelope rotation 




(15) 
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Normal mass loss rate 
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Fig. 5. — Evolution of orbital separation in units of initial orbital separation ao = 5Rg{0) = 
4O5i?0, from the moment synchronization is achieved until end of calculation, as explained 
in the caption of Fig. [Hi Top panel is for the normal mass-loss rate with = 3 x 10~^^ 
in equation ([8]), and the bottom panel is for the high mass-loss rate with rjn = 3 x 10^^'^. 
The initial mass of the RGB star is indicated for each run. All calculations use a companion 
mass of M2 = 0.2Mq, and the RGB core mass at the beginning of the calculation is as in 
the other runs in this section. 
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Fig. 6. — Evolution of the envelope mass for the cases presented in Fig. [5] that enter the CE 
phase. The envelope mass M^nv is scaled to its initial value at the moment of synchronization 

Me„,(0). 
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angular velocity when a system that suffers a Darwin instability just enters the CE phase. 

At the first time that tidal interaction takes place, either a stable synchronization occurs or 
the companion rapidly enters a CE phase. In the case where no stable synchronization occurs, we 
find that when the companion just enters the RGB envelope the envelope angular velocity is 

^ /i /m„ + M2 f Fa^ \ M2 

Ogo = — = — TT — W \ ~B ^ - — ' synchronization, (16) 

OJB ??Mcnv Y Mg \ \ RgO J Menv 

where = MgM2/ {Mg+M2) is the reduced mass, and where in the last equality we have substituted 
our standard values of 77 = 0.2, oq = 5Rg, and took M2 <C Mg. We can substitute Mgnv = 
Mgo — OAMq, and recast equation (fT6]) in the form 

M2 

-^gO — 6— I-O.4M0, no synchronization. (17) 

This is a straight line in the M2 — Mg plane, as evident in Fig. [71 where we present the values of il.go 
when the companion enters the envelope for systems that did not reach stable synchronization. It 
is evident from Fig. [7] that even when the companion spirals-in very rapidly to the RGB envelope, 
it manages to spin-up the primary to high rotation velocity. Such rotation speeds are expected to 
enhance the mass-loss rate, even before considering the gravitational energy that is released by the 
secondary as it spirals-in inside the envelope. 

Consider systems that do reach stable synchronized orbit. The scaled envelope's angular 
velocity when synchronization is first achieved, ilgi, for these systems is presented in Fig. [5J The 
condition to achieve a stable synchronized orbit is given in equation ([7D. (In eq. [7] the assumption 
M2 <C Mg was used, while in the numerical solutions presented here the exact equations are 
used; diff^erences are very small). When such a stable synchronized orbit is achieved, the envelope 
angular velocity, by definition, is equal that of the orbital angular velocity as given by Eq.©. 
It is evident that when synchronization is achieved the primary rotates quite rapidly, and strong 
magnetic activity and enhanced mass-loss rate are expected. Note that more massive secondaries 
bring the RGB envelope to synchronization when they are at larger orbital separations, and hence 
these systems have lower values of il^i. 

As can be seen in Fig. [H most binary systems that have a stable synchronized orbit encounter 
the Darwin instability, and the secondary rapidly spirals-in to form a CE. The normalized RGB 
angular velocity at the moment the system encounters the Darwin instability, ilg2, is presented in 
the top panel of Fig. [9l The orbital separation at that moment is presented in Fig. [TOl As the 
secondary rapidly spirals-in during the instability mode, it transfers orbital angular momentum to 
the envelope, further spinning it up. The normalized RGB angular velocity at the moment the 
secondary of such a system enters the envelope, ^Igs, is presented in the bottom panel of Fig. [9l 
Fig. [9] shows that the primary can reach very high angular velocity prior to the formation of the 
CE, and further strengthens our point that the pre-CE evolution must be considered. 
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Fig. 7. — The envelope angular velocity scaled by the break up velocity when the companion 
enters the envelope iP-go). The systems presented here did not achieve stable synchronized 
orbit at all, and therefore the secondary spiralled in directly into the RGB envelope. The 
black rectangle enclosed the possible area occupied by the progenitor of HD149382 as we 
discussed in section HJ 
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Ugi when synchonization is first achieved for high mass loss rate 




M2(Mq) 

Fig. 8. — The envelope angular velocity scaled by the break up velocity when stable syn- 
chronization is first achieved {^gi)- This figure has the same parameters as Fig. |U but here 
only systems that reach synchronization and form a CE later in the evolution are shown. 
The mass- loss parameter has rjR = 3 x 10~^^ in equation ([H]). 
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flg2 when entering Darwin instability for high mass loss rate 




0.05 0.1 0.15 0.2 

M2(Mq) 

r2g3 when Rg = a for high mass loss rate 




M2{Mq) 

Fig. 9. — The envelope angular velocity scaled by the break up velocity. Top panel repre- 
sents the envelope angular velocity when the Darwin instability occurs {^g2)- Bottom panel 
represents the envelop's angular velocity when the secondary enters the envelope (^2^3). The 
three vertical lines represent the possible locations of the progenitors of the three systems 
discussed in section SJ 
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In should be noted that although the end result of the systems presented in Fig. [7] and in 
Fig. El is the formation of a common envelope, there is a fundamental difference between these 
two evolutionary routes. The systems presented in Fig. [7] form a CE shortly after a strong tidal 
interaction takes place, while those presented in Fig. [9] had evolved through a stable synchronized 
orbit phase. During that phase the RGB star can lose a substantial fraction of its envelope. 

4. DISCUSSION AND SUMMARY 

We have shown that for a large parameters space binary systems reach stable synchronized 
orbit before the onset of a common envelope (CE) phase. Although we have made the calculations 
for RGB stars when their core mass at the beginning of the calculation is Mc{0) = OAMq, the 
results have wider implications. Such stable orbits can be achieved for other evolved stars, like 
AGB stars, and for evolutionary stages earlier or later than those used here. Such stable orbits 
allow the RGB (or AGB, or other evolved giant stars) to lose mass before the onset of the CE. 
Even after the secondary enters the giant envelope, the rotation velocity is high enough to cause 
enhanced mass-loss rate. 

Our results imply that in many cases, studies that neglect the pre-CE evolution (and start 
with a secondary on the surface of a non-rotating giant) might reach some wrong conclusions. For 
example, they might conclude that there is a need for extra energy source to expel the envelope, 
because they do not consider the mass-loss process prior to the onset of the CE phase. 

Our results can be applied to specific systems. We consider here four systems of sdB stars 
with companions: HW Vir with (Mi,M2,a) = {0MMq,0.18Mq,0.89Rq) (Wood et al. 1993); 
PC 1336-018 with (Mi, Ma, a) = (O.5OM0, O.ITMq, O.79i?0) (Kilkenny et al. 1998); HS 0705+6700 
with (Ml, Ms, a) = (O.483M0, 0.134M0, 0.8^0) (Drechsel et al. 2001); and HD149382 (Geier et 
al. 2009) with (Mi,M2,a) = (0.29 - O.53M0, 0.008 - O.O22M0, 5 - 6.1i?0) 

We mark the possible locations of the progenitors of the first three systems on Fig. [9l The 
initial mass of the primary is not well constrained, but these systems most probably went through a 
phase of a stable synchronized orbit. The mass-loss rate could have been high, and a non-negligible 
amount of the RGB envelope has been expelled before the formation of a CE. This increases the 
survivability of the secondary. If a CE was formed only after the RGB envelope expanded, this can 
explain the observations that the core was massive enough to experience a core helium flash after 
the onset of the CE. 

The possible location of the progenitor of IID149382 in the Mg — M2 plane is marked on Fig. 
[71 This system did not reach synchronization, and a CE phase was formed immediately after a 
strong tidal interaction occurred. However, the companion did substantially spin-up the envelope, 
and by that increased the mass-loss rate. This reduced the envelope mass, and by that increased 
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the survival probability of the companion. 
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a , for high mass loss rate 
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Fig. 10. — Orbital separation when the Darwin instabihty occurs (and synchronization is 
lost) for the same systems shown in Fig. [Hi 



